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The zero-temperature phase diagrams of imbalanced two-species Fermi gases are investigated in 
asymmetric optical lattices with arbitrary potential depths, based on the exact spectrum instead of 
the Fermi-Hubbard model. We study the effect of lattice potentials and atomic densities to the fully 
paired Bardeen-Cooper-Schrieffer (BCS) state and particularly the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state. It is found that the increasing lattice potential favors BCS at low densities because 
of the enhanced effective coupling; whereas FFLO is favored at intermediate densities when the 
system undergoes a dimensional crossover. Finally using local density approximation we study the 
evolution of phase profile in the presence of external harmonic traps by merely tuning the lattice 
potentials. 



I. INTRODUCTION 

Ever since the direct observation of phase separa- 
tions of unequal two-component ^Li cold atoms in 
experiments the topic of polarized Fermi gases 

has gained enormous interest from theoreticians in- 
terests in the past few years fs'-'^. Such an imbal- 
anced and attractively interacting fermionic system 
makes it possible to study in the laboratory several ex- 
otic fermionic pairing mechanisms, including Sarma 
or breached pairjjj state and Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state[l2,[I3l with finite pairing mo- 
mentum. However, it is hard for FFLO state to be exper- 
imentally observed in three-dimensional (3D) free space. 
It was proved to survive in a very small region in the 
Bardeen-Cooper-Schrieffer (BCS) side of Feshbach res- 
onance and totally vanish in the BEC side, both for a 
homogeneous Fermi gas@-[l| and a trapped one@. Re- 
cently more and more attention has been paid to one- 
dimensional (ID) [3, [l^ and two-dimensional (2D)(T6j 
systems. It has been pointed out that for these reduced 
dimensions, the polarization window of FFLO state is 
generally broadened compared to 3D free space. 

Despite many previous studies on imbalanced fermions 
in pure ID, 2D and 3D systems, the crossover regimes 
between different dimensions are poorly understood and 
still need further investigations, such as quasi-lD and 
quasi-2D systems. Experimentally, these geometries can 
be achieved by exposing magnetic or optical traps or opti- 
cal lattices along selective directions. Among all trapping 
confinements, optical lattices are the most complicated 
because of the band structures. A very interesting and 
important question then immediately arises: How would 
the dimension and band structure affect the fermionic 
superfluidity, especially the exotic FFLO state? There 
have been several works related to this subject. For ex- 
ample, the FFLO state was studied in pure ID, 2D and 
3D optical latices and shown to be enhanced as a re- 
sult of the Van Hove singularities [l3|; it was also previ- 
ously found with evident proportions in an array of ID 
tubes produced by 2D optical lattices 18[ and in two-lag 



ladders[l9j. In the presence of optical lattices, however, 
all these studies have adopted the Fermi-Hubbard model 
based on the tight-binding approximation (TBA), which 
is generally believed to be valid for deep lattices and 
not very strong interactions such that the lowest band 
model is good enough to describe the system. On the 
other hand, a different approach based on the exact lat- 
tice spectrum ^20] shows that TBA does not necessarily 
produce the correct results because of the multiband ef- 
fects and deviations of the lowest band structures, even 
for optical lattices with considerable band gap. In this 
work we are going to use the exact spectrum to investi- 
gate the zero temperature phase diagrams for imbalanced 
fermions in asymmetric optical lattices, tuned by lattice 
potential depths and atomic densities. This work studies 
the FFLO state in optical lattices with arbitrary poten- 
tial depths. The results obtained reveal the fundamen- 
tal effects of the dimension and the effective coupling to 
the strength of fermionic superfluidities, including FFLO 
state, and thus give an answer to the question previously 
proposed. Therefore our results can serve as a guideline 
for cold fermions experiments and might contribute to 
finally observing the novel FFLO-type pairing state. 

In this work we consider optical lattices applied in se- 
lective directions and take into account the couplings be- 
tween different reciprocal lattice vectors, which rece ntly 
were proved to be necessary for optical latticesfl^ l2ll |. 
The possibility of FFLO pairing is studied in the rema- 
nent free direction(s). We find that by increasing the 
lattice depths, FFLO gradually diminishes at low atomic 
densities as a result of the enhanced effective coupling, 
but revives at intermediate densities where the system 
undergoes a dimensional crossover. Compared with the 
density-driven BCS-BEC crossover initially studied in ex- 
citon superconductors^^ and then applied to attractive 
Fermi gases ^23], the density-driven crossover from weak 
to strong coupling limits mentioned here is unique to op- 
tical lattices and may induce even more significant effects 
to the system in the sense that it changes the effective 
dimension. In contrast, the strength of unpolarized BCS 
superfluidity is enhanced at low atomic densities but frus- 
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trated at intermediate ones because of the discontinuity 
of the energy spectrum. All these results can be un- 
derstood when investigating the density of state (DOS), 
which has an intimate and sensitive dependence on the 
lattice potential and the atomic density. Finally, in the 
presence of the external harmonic trap, we use local den- 
sity approximation (LDA) to investigate the evolution of 
phase profiles merely tuned by the potential depths of 
asymmetric optical lattices. 

The article is organized as follows. In Sec. |TT]we outline 
our general model and method. In Sec. IIIII we give the 
phase diagrams for density-imbalanced two-component 
fermions in several asymmetric optical lattices, including 
quasi-lD and quasi-2D geometries. We mainly focus on 
quasi-lD case, which could be achieved by applying opti- 
cal lattices to either 3D or effective 2D space. The effect 
of external harmonic confinement is studied in Sec. IIVI 
We summarize our results in the last section. 



II. GENERAL MODELS 

We consider a two-component Fermi gas with contact 
interactions by the following Hamiltonian 

g j dYtl;\{Y)il;l{Y)tl;^{Y)il;^{Y), (1) 

where Hq — X]i=a; y z composed by the kinetic en- 
ergy and optical lattice potential if applied; g is the bare 
contact interaction described by s-wave scattering length 
fls in 3D as i = -r^ V„ tt-; and by the binding 

^ g 47ras V 2€q ' ^ 

energy Eb in effective 2D as - = — 4 V„ 75 — tt^- The 
renormalization of g here is to eliminate the unphysical 
divergence due to the high-momentum contribution in a 
Fermi gas. 

To illustrate our scheme, we take a quasi-lD Fermi 
gas confined by two orthogonal optical lattices in x and 
y directions. For example, Hq = J2i=x y z —df/2m + 

[sin^(7ra:;/a)-|-sin^ (Try/a)] , here a is the lattice constant 
and Vq the potential depth. By mean-field treatment, 
first we expand each field operator in terms of eigen- 
wavefunctions of Hq, ipcri^) = Snk <^nk(r)'0nkcr with 

0„,(r) = ^e''=^^^a„4Gx)e*"'-+^-)---, (2) 

here e„k = + e„,fe, -I- Cn^ky (e^^ = 2^) ^^'^ eigen- 
energy; n = {nx,ny} indicates the band index; kj_ lies in 
the first Brillouin Zone (BZ) while has no constraint; 
G± = 2n I a{lx,ly) is the reciprocal lattice vector. The 
Bloch wave functions and eigen-energies can be obtained 
from the decoupled Schrodinger equation in each (x or y) 



direction 

G' i 

^0,nk{G') — ^nkO-nkiG), (3) 

and the eigenvectors satisfy J2g ^nkiG)aii'k{G) ~ (5„„' 
and a„ G) = a^kiG)- For convenience we rescale 
the energy and momentum respectively in units of the 
recoil energy Efi = 2^(f )^ and lattice reciprocal vector 
Then the lattice potential and atomic density are 
expressed by two dimensionless parameters: s — and 

Toward the FFLO state, we employ the simplest single 
plane wave ansatz to the free direction as A(r) = A^e*'^. 
For the lattice part, however, we focus on the most prob- 
able pairing mechanism, i.e., pairing with two opposite 
crystal momenta within the same band. The pairing be- 
tween different bands [l^ is to be neglected here in our 
work. One important reason is that we are dealing with 
the weak coupling regime and therefore the interaction is 
not large enough to overcome the energy differences and 
form inter-band pairs. 

By employing two set of pairing fields in terms of Bloch 
state indices (nk) and reciprocal lattice vectors (Q), 

nk 

A„k - E^Q^^nk*> (4) 

Q 

with A/^i^ = EG«n-k(-G)a„k(G + Q) and G, Q all 
lying in x-y plane, the Hamiltonian is then reduced to 

cr nkcr 
nk Q 

By diagonalizing the Hamiltonian, the thermodynamic 
potential is calculated at zero temperature as 

^ = :^5l{®(-^"k+)i^nk++e(-£;nk-)£;nk- + 
nk 

Ck-VCk^ + A„,^}-E^' (5) 

Q ^ 

where the quasi-particle spectrum reads 

^^nk±- VCk'+A„k'±Ck (6) 
with Ck = ^"xfc. + ^nyky + i^t+qj2 + e-fe.+q./2)/2 - M> 

Ck = i^t+q./2 - '=-fe.+g./2)/2 - h, and chemical poten- 
tials /i = (;zt + Ai4,)/2, {fif- fii)/2. 
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From the saddle-point equations cJrj/SAq = 0, 
dn/dqz = and N„ = —dQ/dfi^r, we obtain the gap, 
current and density equations as 



- E 



'nk ~ nk 



nq.l1 + E - E 



(7) 
(8) 



B„k+<0 B„k-<0 



= ^( E 1- E !)■ (9) 



S„k+<0 B„k-<0 



Therefore the consideration of couplings between dif- 
ferent Q directly results in the coupled gap equations 
d?]). To get solutions with great precision one must take 
into account as many non-zero Q as possible, and we find 
that the smaller ones dominate over the larger ones, es- 
pecially in free space one has Aq = Aq^qo- This could 
be understood when examining the properties of M^^,, 
which is closely related to the pairing amplitude Aq for 
each Q, for a simple ID optical lattice. At very high 



energy levels (n » 1), M, 



Q 



JO close to those in free 



space, so the contributions mainly come from several low 
bands. Near the bottom of the lowest band, the pertur- 
bation theory justified for shallow lattices (s ^ 1) gives 

^nfc = 1(Q = 0), |(Q = ±1), ^(g = ±2) and zero 
for other high-order ones; whereas the TBA for deep lat- 
tices (s 3> 1) gives a smooth Gaussian distribution as 
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Q ). For intermediate s it is numeri- 
cally verified that smaller Q always lead to predominantly 
larger M^^, and therefore contribute the most to the sum- 
mations in Eq.(|3]) and ©. 

In view of the properties of non-zero Q— pairing in 
this case, besides Q = we consider the other four 
smallest nonzero ones: (±1, 0), (0, ±1). Because of the 
equivalence of x and y directions, all these non-zero Q 
share the same pairing amplitude Ai. Therefore we get 



Ank = Ao + 2Ai(A4Jf ^* + A/^k ^*)' and the coupled gap 
equations are in terms of Aq and Ai. 

Next we utilize the this model to study the phase dia- 
grams of imbalanced fermions, either in terms of chemical 
potentials {h,fi) or in terms of particle densities {Sn,n). 
In the /i— /i phase diagram, the phase boundaries separat- 
ing the fully paired BCS state, FFLO state, and normal 
(N) state are determined as follows. We solve Eq.Q for 
BCS state at given /i and obtain the minimum excita- 
tion energy Emtn = niin y^Clfe^ + ^nk the first step. 
Then the BCS-normal (BCS-N) phase boundary {Hen) 
is determined from the following equation 



(01)* 



^BCsilJ-. hBN, A = {Ao, All) = hBN, A = {0}). 

(10) 



Note that this equation has a solution for hBN{< Emin) 
only when ilscsih = Emin) > ^Nih = Ej^in)- Oth- 
erwise the stable magnetized superfluid {SFm) phase 
would interpolate between BCS and normal together 
with two continuous phase boundaries in between. Near 
the FFLO-N boundary, 57 can be expanded in terms of 
small pairing amplitudes as = X]q q' ^Q Q' Aq Aq/ 
where 
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V 



^'^nk^'^nk 



(11) 



For particular pairing momentum (q), the upper limit of 
h for FFLO state (he) is determined by setting \W\ — 0, 
and in this particular case 1^ is a 5 x 5 matrix. While 
the FFLO-N boundaries {hpN and qpN) are obtained by 
finding the maximum of he as a function of q. Under 
mean-field treatment, the BCS-N phase transition is al- 
ways found to be first-order, while FFLO-N is second- 
order. Moreover, as the thermodynamic potentials of 
FFLO and normal states are very close to each other 
near the BCS-FFLO boundary {hsp), which makes hsp 
very close to h^Ni in the following we do not distinguish 
these two boundaries. 

After obtain the h — fi phase diagram, we directly con- 
vert it to the polarization-density (P — n) phase dia- 
gram composed of phase separation (PS) of BCS and 
N state, FFLO and N state. For example, the BCS-N 
phase boundary of the former (Hen) corresponds to the 
PS-N boundary of the latter (Ppat); in other words, PpN 
is just the polarization of a N state with chemical poten- 
tials (/i, ft-BAf), for the reason that the PS with maximum 
polarization just represents the critical point when BCS 
state vanishes and N spreads to the whole space. Simi- 
larly the FFLO-N boundary of the former (/ifw) corre- 
sponds to the same one of the latter {Pfn)- 



III. PHASE DIAGRAMS 

The phase diagrams of imbalanced two-component 
fermions including FFLO state have been studied previ- 
ously in free 3D 0431 and 2 p[l6ll spaces based on mean- 
field theory, and in a 1 D [T5|space using Bethe ansatz 
technique. Remarkably in ID case FFLO exhibits a quite 
broad polarization window ranging from a rather small 
value to unity in the weak coupling limit. In this section 
we study the quasi-lD and quasi-2D geometries that can 
be achieved by applying optical lattices in selective di- 
rections. Note that there are two ways to generate the 
quasi- ID system. One is from the combination of a tight 
harmonic confinement in one direction and an optical lat- 
tice in another, which generate ID tubes lying in a plane; 
the other is from two orthogonal optical lattices, which 
induce ID tubes in 2D lattice sites. In the next sub- 
sections we illustrate these two cases. Finally we briefly 
introduce the quasi-2D system generated by a ID optical 
lattice. 
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For realistic simulations, we apply the cutoff momen- 
tum well above the Fermi momentum kp, which is de- 
fined for an unpolarized normal Fermi gas with the same 
total atoms number as the polarized one studied, to en- 
sure accuracy. The attractive s-wave interaction is fixed 
to be well within the weak coupling limit. 



A. Quasi-ID geometry in an effective 2D space 

For an effective 2D space with the axial freedom of 
motion frozen by a tight harmonic confinement, the 
renormalized atom-atom interaction is characterize d by 
the two-body binding energy Eb = ^exD(\/27r j-) [25|, 

where w is the confinement frequency, I — \J^^^ is the 
characteristic confinement length and C ~ 0.915. In such 
a 2D systemj_all phase boundaries could be analytically 
obtained[l^ except for FFLO. We numerically cal- 
culate the FFLO-N boundary and verify the previous 
predictions [13, [11] in weak coupling limit that two Fermi 
surfaces just touch at the critical point, with a constant 
pairing momentum amplitude q ~ 2y/mEB- The two 
relevant boundaries are 



hBN - 

^EB{2fi-EB). ^i>Eb 



A* > — 2 — 



(12) 
(13) 



The phase diagrams for a 2D Fermi gas are shown 
in Fig. [TJ and here three points are emphasized in or- 
der. First, numerical simulations show that the FFLO-N 
boundary hpN, which is the maximum value of he at cer- 
tain non-zero q, shows up immediately when /i > 0, but 
is less than the BCS-N boundary Hbn until /i up to 
at point C [shown in Fig. Hl^a)]. Only beyond C, FFLO 
can exist as a candidate of ground state. Second, Fig. 
[ijb) shows that the maximum polarization window for 
FFLO state to exist as ground state is less than 20%, even 
when taking into account the small distinction between 
hsF and Hbn- Third, throughout the range of /i, there 
is no stable magnetized supcrfluid (SFm) phase found as 
that in the strong coupling limit in 3D free space. One 
powerful piece of evidence is that the Sarma-N boundary 
(hsN) can never exceed the lowest BCS excitation energy 



hs 



N 



Er, 



E b/2, < Eb/2 

Vf^EB/2, fi>EB/2 

/i + Eb, -Eb /2 <fi<0 
^EBi2^i + EB), M>0 



(14) 
(15) 



implying the exclusion of SFm from the ground-state 
solutions. Here hsN is obtained by solving the gap equa- 
tions at zero gap amplitudes. We also use this criterion to 
judge the existence of SFm in the following calculations 
for optical lattices. 
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FIG. 1: (color online) Phase diagrams of a 2D Fermi gas. 
(a)/i — /i diagram constructed by vacuum, BCS, FFLO, par- 
tially polarized normal (PPN) and fully polarized normal 
(FPN) states. The boundary hsF is just slightly below Hen- 
Red cross points Esyh/ Eb) from left to right: ^(1,1), 
B(^^,^^) and C(|, ^f) correspond to hpN = A*, 
hsN = l-i. and hpN = fesjv respectively. Orange dashed line is 
for h = (b)P — n diagram constructed by phase separa- 
tion (PS), FFLO and normal (N) states. Inset shows two AP 
plots APi = Pfn — PpN and AP2 = Pfn — Ppf, with each 
peak {fj./EB, APmax) locating respectively at _D(3. 605, 0.146) 
and S(3.750, 0.177). 



When applying an optical lattice along one direction 
the phase diagrams are dramatically changed as shown 
in Fig. [21 First we study the low density limit. It is 
observed that the minimum density for FFLO to exist, 
rimin, moves to the right (becomes large) as lattice po- 
tential s increases. We argue that this indicates a much 
more enhanced coupling of atoms moving in deeper lat- 
tices. Actually this critical point {rimin) is an analog 
of point C in Fig. [Ija) for free 2D system, in which 
nmin monotonously increases with ^ic{= 5Eb/4:) and 
also Eb- Optical lattices applied here drive the system 
to a strongly interacting regime characterized by the en- 
hanced effective binding energies and thus increase n„iin- 
Moreover, in this dilute limit the maximum of critical 
chemical potential differences hc{q) is always at q = 
(see Fig. Efa)), indicating no FFLO but only an unsta- 
ble Sarma state {hsN < Emm)- 

Second we observe that as the density n or /i increases, 
the FFLO-N boundary hpN initially goes up but then 
suddenly drops at certain position, the reason for which 
is illustrated in Fig. [31 As n increases, hpN with nonzero 
pairing momentum q starts to exceed Hbn which stabi- 
lizes FFLO states. Meanwhile, the shape of Fermi sur- 
face distorts from an ellipse to two disconnected lines, 
implying a crossover from 2D to a quasi-lD geometry 
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FIG. 2: (Color online) Phase diagrams of a 2D Fermi gas 
when applying optical lattices along an arbitrary direction 
with different depths s = 3, 5. The gas is in the weak inter- 
acting regime with the two-body binding energy Eb = 0.2_E_r. 
(a)/i — Afi diagram with A/i calculated from the bottom of 
the lowest band. Orange dashed line is for h — |A/i| defining 
the threshold for a fully polarized normal state. BCS, FFLO 
and normal states are separated by the corresponding bound- 
aries as labeled in the legend. FFLO exists as the ground 
state when Hen < h < hpN- (b)P — n diagram. From the 
bottom to top are phase seperation, FFLO and normal state 
in turn. At a certain density, Pfn suddenly drops to a lower 
value corresponding to the same behavior of hpN in (a). 



and correspondingly hp^ has a rapid increase during the 
crossover. When n increases further, the shape of Fermi 
surface changes again when atoms begin to fill a higher 
lattice band. However, this occupation of a higher band 
exposes a destructive effect on FFLO pairing along the 
free direction. For example Fig. shows that the max- 
imum value of he for Afi/Eji = 2.73 is much less than 
that for 2.43, leading to a sudden drop oi hpN as well as 
PpN in Fig. [2] Moreover, it is noticed that the discon- 
tinuity of hpN in this sensitive region produces double 
values of Pfn at a definite density n as a side effect, 
in which case the larger one should be chosen as the real 
critical polarization for FFLO. Finally it is shown by Fig. 
Wlb) that at intermediate n the PS-N boundary PpN is 
suppressed by increasing s, in contrast to Pfn- We at- 
tribute this to the discontinuous lattice spectrum that 
reduces the Hilbert space for pairing and therefore disfa- 
vors the BCS superfiuidity. 

Finally we comment that these effects of optical lattices 
at low and intermediate densities can be well understood 
from the point of view of density of state (DOS), see 
also Fig. m The variations of DOS bring two opposite 
effects to the BCS superfiuidity (SF) and FFLO. At low 
densities (n), DOS is enhanced by the optical lattices 
applied and therefore favors SF but disfavors FFLO as 



FIG. 3: (Color online) (a)Critical chemical potential differ- 
ences (he) as a function of total pairing momenta (q) at sev- 
eral averaged chemical potentials (A/x) with lattice potential 
s = 5. (b)Fermi surface contour in momentum space. The op- 
tical lattice is applied along x direction. The inner and outer 
red dashed lines respectively represent the Fermi surfaces at 
the top of the lowest band (labeled by 0.264) and slightly be- 
yond the bottom of a higher band (2.73). The outer line also 
corresponds to the dashed-dotted curve in (a). 
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FIG. 4: (Color online) Density of state (DOS) scaled by 
1/(0^ Er) at the Fermi surface versus particle densities n = 
No? /S. Optical lattices are applied along an arbitrary di- 
rection in 2D space with different depths s = 0, 3, 5. In free 
space {s = 0), DOS is a constant (7r/4). 



a result of the strong effective couplings; at intermediate 
n, DOS decreases with n which approximately satisfies 
DOS ~ 1/n indicating the evolution to ID geometry, so 
this time FFLO is favored while SF is suppressed. When 
n further increases, DOS becomes stable at the same level 
as in 2D free gas when atoms start to fill a higher lattice 
band. After this point the quasi-lD geometry fades away 
and FFLO is no longer favored. 
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B. Quasi-ID and quasi-2D geometries in a 3D 
space 

The quasi- ID geometry in 3D space is formed by ap- 
plying optical lattices along two orthogonal directions. 
The theoretical model has been outlined in Sec. HIl and 
here we present the phase diagrams in Fig. [Sj 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
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FIG. 5: (color online) Phase diagrams of a 3D Fermi gas when 
applying optical lattices along two orthogonal directions with 
different depths s — 3,5. The s-wave scattering length is fixed 
to be as = —a/3. {a.)h — A/i diagram. Orange dashed line is 
for h = |A/^|. (b)P — n diagram. Each boundary separates 
corresponding phases in the same way as in Fig. [H 

At a given low density (n), the maximum polarization 
for phase separation (Ppn) increases from a small value 
to unity with lattice potential s, in contrast to quasi-lD 
geometry in 2D space; however it does reveal the same 
fact that optical lattices enhance the effective coupling. 
Moreover in this region, hp^ < hpN so actually FFLO 
will not serve as a candidate for ground states. As n con- 
tinues to increase, the system undergoes a crossover to 
quasi-lD geometry, and hp]\j gradually goes up. Com- 
pared with a negligible FFLO polarization window for 
the 3D free Fermi gas, the optical lattice with s ~ 5 can 
broaden the window from as large as 25% at the optimum 
density n ~ 2.5. When n or /i is large enough to touch 
a higher lattice band, hpj^ and PpN again suffer from a 
drop behavior as introduced in the previous section. 

Next we turn to quasi-2D geometry in 3D space which 
could be generated by an axial optical lattice. We do not 
present the phase diagrams here but give a brief introduc- 
tion. Compared with quasi-ID geometry the available 
FFLO window in this case is much less robust, for in- 
stance less than 10% with lattice potential s — 5. More- 
over there is no obvious sudden drop of hpN or PpN. 
This might be ascribed to two tunable parameters {\q\ 
and 9) in this case for FFLO pairing momentum vector 



{<lxi<ly) = \q\{cos9,sm9), which make the pairing rela- 
tively flexible and easily adaptive to external variations. 



IV. PHASE PROFILE UNDER EXTERNAL 
HARMONIC CONFINEMENTS 

In this section, we study the effect of asymmetric lat- 
tices to the phase profile of imbalanced fermions in a 
harmonic trap using local density approximation (LDA). 
In LDA, locally the gas is considered to have the same 
properties as a bulk gas, which depend on the local av- 
eraged chemical potentials /Lt(r) = (/.tot + Mo^)/2 — V^(r) 
and position-independent difference h = {^q^ — ^oi)/2. 
MO(t(o' =t,4-) is the chemical potential of spin-cr atoms at 
the trap center, and can be self-consistently determined 
by the total particle numbers N^, s-wave interactions and 
lattice potentials. 




FIG. 6: (color online) . (a, b,c) Spatial distributions of the nor- 
malized total density (n/n(0)), density difference (dn/n(0)), 
gap amplitudes (Ao,Ai), FFLO pairing momentum for 
a density-imbalanced two-species ^Li Fermi gas. n(0) = 
N{0)a^ /V is the dimensionless total density at the trap cen- 
ter. The optical lattices s = 0(a), 3(6), 5(c) are applied along 
X and y directions with lattice constant a = 532nm. An 
isotropic harmonic trap is also applied with trapping fre- 
quency / = 200Hz. The s-wave scattering length is fixed 
to be tta = —a/3. Other parameters are chosen so that 
these three plots (a,b,c) have nearly identical particle num- 
bers: s,N/2,P,h/ER,n{0) = 0, 4.3(10**), 0.42, 0.31, 3.03(a); 
3, 4.35(10^), 0.42, 0.32, 2.44(6); 5, 4.45(10^), 0.43, 0.34, 2.35(c). 
(d) Spatial distributions of the averaged total density (n/2) 
and polarization (P) corresponding to (a-c). 

In Fig. [S]we give three typical phase profiles under an 
external isotropic harmonic trap. The optical lattices are 
applied in two orthogonal directions and form an array 
of ID tubes in 2D lattice sites. Using LDA, we solve the 
ground state at each position by minimizing the thermo- 
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dynamic potential ri(Ao, Ai, Qz) in terms of its three pa- 
rameters. Meanwhile /igt ^'^d are adjusted such that 
the particle numbers are almost identical for different s. 
Therefore Fig. [5] actually shows how the phase profile of 
an imbalanced Fermi gas in a harmonic trap evolves when 
switching on the optical lattices, especially for the FFLO 
state. Without lattices (s = 0), the profile is mostly con- 
structed by BCS and normal state, with negligible FFLO 
window in between. While increasing s, the window is 
gradually broadened. For s = 3, the spatial range of 
FFLO has extended to be nearly as the same as that 
of BCS and normal state. During this period, particles 
move from center and edge to the middle region, causing 
both n(0) (density in the trap center) and Rtf (Thomas- 
Fermi radius) to decrease correspondingly. By increasing 
s further, FFLO would gradually spread to the trap cen- 
ter and take over BCS state, until finally only FFLO and 
normal state are left in the trap, see Fig. [B^c) for s — 5. 
Moreover, notice that the discontinuities of densities ng., 
gap amplitudes Aq, Ai at the BCS-FFLO boundary im- 
ply a first-order phase transition, while the continuities 
at FFLO-N boundary imply a continuous transition. 

Finally two statements about FFLO are given as fol- 
lows. First, FFLO should be observed far away from the 
trap edge due to the atomic dilution and strong effective 
interactions there; second, its spatial range is enlarged 
by the increasing geometric asymmetry. 



V. SUMMARY 



In conclusion, we have studied the phase diagrams of 
imbalanced two-species fermions in asymmetric optical 
lattices. We found that the optical lattices applied ex- 
pose two opposite effects on the BCS and FFLO states, 
depending on the atomic densities. One is to enhance the 
effective interaction at low densities which favors BCS 
but destroys FFLO; the other is to induce the crossover 
to lower dimension at intermediate fillings which favors 
FFLO, but destroys BCS due to the discontinuity of lat- 
tice energy spectrum. Among all the asymmetric geome- 
tries, we find that the quasi- ID system is the most fa- 
vorable one for observing the FFLO state. Finally by 
using LDA, we present several typical phase profiles in 
a harmonic trap with substantially different FFLO pro- 
portions, which are merely tuned by the potential depths 
of applied asymmetric optical lattices. These lattice ef- 
fects still need to be further explored in cold-atom exper- 
iments. 
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